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ABSTRACT 

Accurate maps of Galactic reddening are important for a number of applications, such as mapping the 
peculiar velocity field in the nearby Universe. Of particular concern are systematic errors which vary 
slowly as a function of position on the sky, as these would induce spurious bulk flow. 

We have compared the reddenings of Burstein & Heiles (BH) and those of Schlegel, Finkbeiner & 
Davis (SFD) to independent estimates of the reddening, for Galactic latitudes \b\ > 10°. Our primary 
source of Galactic reddening estimates comes from comparing the difference between the observed B — V 
colors of early-type galaxies, and the predicted B — V color determined from the B — V- Mg2 relation. 
We have fitted a dipole to the residuals in order to look for large-scale systematic deviations. 

There is marginal evidence for a dipolar residual in the comparison between the SFD maps and the 
observed early- type galaxy reddenings. If this is due to an error in the SFD maps, then it can be corrected 
with a small (13%) multiplicative dipole term. We argue, however, that this difference is more likely to be 
due to a small (0.01 mag.) systematic error in the measured B — V colors of the early-type galaxies. This 
interpretation is supported by a smaller, independent data set (globular cluster and RR Lyrae stars), 
which yields a result inconsistent with the early-type galaxy residual dipole. BH reddenings are found 
to have no significant systematic residuals, apart from the known problem in the region 230° < I < 310°, 
-20° < b < 20°. 

Subject headings: dust, extinction — galaxies: distances and redshifts — galaxies: elliptical and 
lenticular, cD — globular clusters: general — stars: variables: other 



1. INTRODUCTION 

Accurate values of Galactic extinction are essential for 
a number of applications, such as measuring distances and 
peculiar velocities of galaxies. A small error of only 0.015 
mag. in the R-band Galactic extinction corresponds to a 
systematic distance error of 1% for the Fundamental Plane 
distance indicator. When measuring the large-scale bulk 
motion of galaxies, the small-scale accuracy of reddening 
maps is not critical, since random errors add to the scatter 
but do not introduce a systematic bias in the bulk flow. On 
the other hand, large-scale coherent systematic errors in 
the reddenings would introduce corresponding systematic 
errors in the bulk flow. 

In this paper, we test two reddenings maps for such 
large-scale systematic errors: the reddening maps of 
Burstein & Heiles (1982, hereafter BH), which are based 
on neutral hydrogen (for Dec. > —23, Shane- Wirtanen 
galaxy counts are used to determine dust-to-gas ratios); 
and the maps of Schlegel, Finkbeiner & Davis (1998, here- 
after SFD) which are based on dust emission measured 
by IRAS with dust temperatures determined from DIRBE 
data. Note that the two maps are completely independent. 

There is some evidence that the BH reddenings are sys- 
tematically in error in at least one large coherent region of 
the sky. Burstein et al. (1987, hereafter B87) found that 
the BH reddenings were overestimates (by a factor 2) in 
the region 230° < I < 310°, -20° < b < 20° (hereafter 
referred to as the Vela region). Note that this region lacks 
a dust-to-gas ratio estimate. 

SFD calibrated their reddening maps using the colors of 
brightest cluster galaxies from Postman & Lauer (1995) 



and tested their maps with the B — V- Mg2 relation and 
the data of Faber et al. (1989). However, SFD did not 
test for coherent residuals as a function of longitude. Such 
residuals are the bane of peculiar velocity work; they are 
the focus of this paper. 

To test the BH and SFD maps, we require a sample of 
distant objects with standard colors. Our primary sample 
is the early-type galaxy data set of Faber et al. (1989). 
We will make use of the tight correlation between Mg2 
index and B — V color. Lynden-Bell et al. (1988) used 
the B — V- Mg2 relation to show that putative systematic 
errors in the BH reddenings could not be large enough 
to explain the bulk flow of 570 km s _1 towards I = 307°, 
6 = 9°. However, Lynden-Bell et al. did not attempt to 
push this farther and obtain constraints on the accuracy 
of the BH reddenings. This was perhaps due to systematic 
uncertainties in the Mg2 index, which are at the level of 
~ 0.01 in Mg 2 (Davies et al. 1987), equivalent to 0.008 in 
E(B-V) or 0.03 in B-band extinction. 

We have recently determined the corrections required to 
bring Mg2 indices from 23 different spectroscopic systems 
onto a standard system, by comparing 1784 repeat Mg2 
measurements for 418 galaxies (Smith et al. 1997, Hudson 
et al. 1999). As a result the uncertainties in the correc- 
tions to the 7S Mg 2 systems have been reduced by a factor 
5. This allows us to re-examine the reddening maps with 
increased accuracy. 

There has been much discussion in the literature of 
the value of reddening at the poles, i.e. the zero-point 
of the reddening maps. This zero-point is impossible to 
calibrate using only extragalactic color standards. Fortu- 
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nately global zero-point errors have no effect on peculiar 
velocity work since they are absorbed into the zero-point 
of the distance indicator. We make no attempt in this pa- 
per to determine the correct zero-point of the reddening 
maps. 

The outline of this paper is as follows. In sj2| we use 
early-type galaxies as redd ening esti mat ors, and present 
our data and method (§^J), results (§ |2.2| ) and disc uss the 
effect of photometric errors in the B-V data (§|J). In §| 
we constrain the systematic errors in the SFD maps using 
independent reddening estimates obtained from globular 
clusters and RR Lyrae stars, and we summarize our re- 
sults in §[| 

2. REDDENING FROM EARLY-TYPE GALAXIES 

2.1. Data and Method 

Our primary data are the early-type galaxies of the 
'Seven Samurai' (7S) survey (Faber et al. 1989), with 
B — V colors measured within 67 arcsec from B87. How- 
ever, we use revised Mg2 data for these galaxies which are 
corrected for offsets between Mg2 measurements obtained 
from different telescopes (Hudson et al. 1999). Our cor- 
rected Mg2 values agree well with the subset of the 7S 
data recently published by Trager et al. (1998). We use 
only the B — V data with D n photometric quality of 1, 
corresponding to random and systematic errors less than 
0.02 dex (B87). The data are given in Table @. Only 
the first and last four lines of the data are reproduced 
in the printed version. The full table is accessible elec- 
tronically. Finally, we limit our analyses to Galactic lat- 
itudes \b\ > 10. The B—V- Mg2 relation has a scat- 
ter of only 0.028 mag. This relation is shown in Fig. [l| 
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Fig. 1. — The B — V- Mg2 relation with B — V corrected using 
SFD reddenings. Four galaxies with Mg2 < 0.15 are not shown. 

If we compare the B — V- Mg2 relations of the subsample 
of cluster galaxies to the subsample of galaxies in groups 



(with velocity dispersions [Hudson 1993] less than ~ 400 
km s _1 ) or in the field, we find that the cluster galaxies are 
redder than group/field galaxies by 0.01 ± 0.004 mag at a 
given Mg2. Hudson (1999) shows that Malmquist bias ef- 
fects can cause offsets between field and cluster samples 
at this level, even when distance-independent measures 
(B — V and Mg 2 ) are compared. For the purposes of this 
paper, however, we are not concerned with the interpreta- 
tion of this offset. We simply correct the B — V colors of 
field galaxies by adding 0.01. The conclusions of this pa- 
per, however, are not affected if we neglect this correction. 

The residuals from the B — V- Mg2 relation for both 
SFD and BH reddenings are plotted on the sky in an Aitoff 
projection in Fig. |[ No obvious trend to the residuals is 
immediately apparent although the Vela region stands out 
in the BH plot, and to a lesser extent in the SFD plot. 
Clearly any deviation of the residuals from isotropy is at 
a fairly low level. To measure the anisotropy requires a 
specific model. 

Our method is straightforward: we use the tight intrinsic 
relation to predict the true B — V and compare this with 
the observed B — V h s , to obtain a reddening estimate 



E(B-V) Emp = (B-V) ohs - [A + B Mg 2 ] 



(1) 



where A and B are free parameters. Since our goal is to 
look for coherent residuals from the SFD and BH maps, 
we fit the map reddening to these reddening estimates al- 
lowing for a correction term. For the form of this correc- 
tion, we adopt the lowest order spherical harmonic, i.e. 
a dipole. We have considered both additive and multi- 
plicative dipole corrections to the reddenings. We adopt 
a multiplicative correction, which is perhaps more realis- 
tic (e.g. a coherent error in the dust temperature). The 
reddening estimates are then fit by 



E(B -V)jwp = {l+F-r) E(B-V) 



map 



(2) 



where E(B — V) map is the predicted reddening given by 
BH or SFD. F is a dipole vector and f is the unit vec- 
tor towards the galaxy. In practice, these two steps are 
performed simultaneously. 

The residuals around these relations are slightly non- 
Gaussian, with a tail of blue objects. In order that our 
fit be robust, we minimize the sum of absolute deviations 
of the residuals. This is equivalent to fitting the median 
(Press et al. 1992). We estimate the significance of the 
dipole coefficients by performing Monte Carlo simulations 
in which the positions of early-type galaxies are random- 
ized while their B — V and Mg2 values are kept fixed. 

2.2. Results 

The result of the dipole fits to both SFD and BH 
are given in Table g. For the SFD Mg 2 fit, we de- 
tect a marginally significant dipole (at the 92% CL) of 
|F SFD | = 0.13^0 towards I = 80°, b = 23° (errors corre- 
spond to 66% confidence range). Note that the anti-apex 
of this dipole is I = 240°, b = -23°, which is in the Vela 
region. 

We detect no significant residual dipole in the BH red- 
denings, if the Vela region is excluded from the analy- 
sis. (If we exclude this region from the SFD analysis, the 
significance of the SFD dipole term drops to 87% CL). 
The lower significance of the BH dipole compared to SFD 
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Table 1 
Early- type galaxies 



Name 


B- 


-V 




1 




b 


Mg 2 


Qual 


E(B-V) SFD 


E(B-V) BH 


Clus a 


Vol* 

km s 1 


A1656:D-107 





.87 


56. 


.7953 


88 


.0890 


0.2311 


111 


0.0118 


0.0114 


1 


821 


A1656:D-125 





.92 


56. 


.9377 


87 


.9992 


0.2531 


111 


0.0089 


0.0121 


1 


821 


A1656:D-136 





.90 


60. 


.1584 


88 


.1393 


0.2672 


111 


0.0105 


0.0120 


1 


821 


A1656:D-153 





.95 


58. 


,7005 


87 


.9678 


0.2747 


111 


0.0079 


0.0124 


1 


821 


N7626 


1 


.04 


87. 


.8581 


-48 


.3781 


0.3260 


111 


0.0725 


0.0411 


1 


212 


N7768 





.94 


106. 


.7153 


-33 


.8139 


0.3153 


212 


0.0386 


0.0276 


1 


-999 


N7785 


1 


.00 


98. 


.5380 


-54 


.2775 


0.2915 


212 


0.0564 


0.0420 








U03696 


1 


.01 


168. 


,5668 


22 


.8446 


0.3011 


312 


0.0715 


0.0716 


1 


327 



"Flag for cluster/group (set to 1) or field (0) 
6 Cluster/group velocity dispersion (—999 = no data) 




Fig. 2. — Early-type galaxy residuals E(B — V) e iiip - E(B — V) m ap are shown on the sky in Aitoff projection in Galactic coordinates. Open 
circles show galaxies with positive residuals, i.e. galaxies for which the maps underprcdict the reddening; crosses show galaxies with negative 
residuals. The size of the symbol scales with the absolute value of the residual. The symbols in the lower left corner show the sizes of symbols 
corresponding to residuals of 0.06 and 0.03 mag. The top panel shows the residuals from the SFD maps; the poles of the SFD correction are 
marked by the plus and minus symbols. Celestial North and South poles are indicated by 'N' and 'S'. The solid line is the equator. The lower 
panel shows the BH residuals. The dashed line is at Dec. = —23° and divides the region with galaxy counts from the region without. The 
dashed box shows the Vela region. 
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Table 2 
dlpoles of e(b — v) residuals 



Sample 




scatter 


Map 


Corr a 


W\ 


I 


b 


Prob : 


Ellip 


311 


0.028 


SFD 


0.000 


0.13 


80 


23 


8.4 


Ellip c 


296 


0.031 


BH 


0.000 


0.17 


64 


17 


44.7 


Ellip c 


296 


0.028 


SFD 


0.000 


0.13 


94 


15 


13.1 


Ellip 


311 


0.028 


SFD 


0.010 


0.07 


328 


-69 


57.2 


Ellip c 


296 


0.031 


BH 


0.010 


0.11 


47 


-27 


87.7 


GC+RR 


137 


0.028 


SFD 


0.000 


0.08 


312 


-47 


95.0 


Ellip+GC+RR 


448 


0.028 


SFD 


0.000 


0.05 


98 


-35 


38.6 



17 



"Correction applied to Northern E(B — V) data 
^Probability that measured dipole could arise by chance 
c Excluding the Vela region 



should not be misinterpreted as evidence that its residual 
dipole is lower - it is not - rather the reduced significance is 
due to the larger scatter in the residuals when BH redden- 
ings are used, which leads to a larger error in the dipole: 
|F BH | = 0.17t°#. 

Fig. |3| shows the residuals from the B — V- Mg2 rela- 
tion as a function of the cosine of the angle from the apex 
(I = 80°, b = 23° for SFD and I = 64°, b = 17° for BH). 

2.3. Discussion 

We have modeled the residuals with a dipole term. How- 
ever, it is also possible that the residuals are clustered 
in a few, smaller-scale (30-60°) patches on the sky. We 
have already discussed the Vela region, in which both the 
BH reddenings as well as those of SFD overestimate the 
reddening. The region 120° < I < 180°, 0° < b < 40°, 
(hereafter referred to as the Cameleopardis region) also 
stands out in Fig. |^ due to its excess of red galaxies. The 
Cameleopardis region is opposite on the sky to Vela re- 
gion, so it is possible that these two regions together are 
the source of the dipole signal. The sparseness of the data 
do not allow us to determine whether the residuals are con- 
fined to smaller areas, or whether they vary more slowly 
(like a dipole). 

At low Galactic latitudes, systematic errors are expected 
to be largest, both in the predictions and in the observed 
colors, due to the problems of contamination by fore- 
ground stars in the Milky Way. If we limit the analysis 
to Galactic latitudes \b\ > 20°, however, we still find an 
identical dipole for SFD and a somewhat larger dipole for 
BH, compared to the results found above for |6| > 10°. 

So far we have interpreted these residuals as errors in 
the reddenings. It is possible that systematic errors in the 
Mg 2 data or B — V data are responsible for these residuals. 
The Mg2 data for the 7 spectroscopic systems tabulated 
by 7S were corrected for systematic offsets, but the uncer- 
tainty in these corrections was 0.01 mag. (Davies et al. 
1987). By inter-comparison of common data, these Mg2 
data have been re-corrected and brought onto a common 
system (Hudson et al. 1999). While the corrections are 
large (as much as 0.02), the uncertainties in these correc- 
tions are now considerably smaller (typically 0.0025 mag). 
This translates to only 0.002 mag in E(B — V) which is 

1 During the K3 run a can of oil was found after the end of the run 



a factor of 5 smaller than the residual dipole. Therefore 
systematic errors in Mg2 cannot account for the residuals. 

It is also possible that there are systematic errors in the 
B — V colors. We have performed fits in which we allow for 
a systematic error in color between galaxies in the north- 
ern and southern celestial hemispheres. In Table ||we also 
show fits for BH and SFD after a 0.01 mag correction ap- 
plied to the Northern B — V values. When this correction 
is applied, the residual dipole terms are small and consis- 
tent with zero. Thus an offset in photometry can account 
for the residuals in the SFD comparison. In fact, the sig- 
nificance of this photometric correction is somewhat higher 
than the dipole term fitted above. 

As discussed above, we cannot really distinguish be- 
tween photometric errors in small patches of the sky and 
large-scale (north-south) photometric offsets. The B — V 
colors tabulated in B87 arise from 10 different photomet- 
ric systems. Table 4B of B87 indicates that, while most 
of their photometric systems are of high accuracy (bet- 
ter than 0.01 mag. in B — V), there are several excep- 
tions. Most notable of these is the K3 system (KPNO 
0.9m Oct 81) 1 . The K3 system shows a systematic offset 
of 0.02 mag in B — V and possibly a weak (uncorrected) 
color term (B87). Note that about half of the galaxies in 
the Cameleopardis region have B — V colors from the K3 
system, so this may account for some of the deviant galax- 
ies in this region. It is beyond the scope of this paper to 
extract the K3 observations from the merged photometry 
of B87. We simply note that although one can interpret 
the residuals as a north-south photometric offset, it is also 
possible that errors from one or more subsystems could 
mimic the same effect. 

Note that in Fig. || both BH and SFD have residuals 
in the same locations (e.g. in the Vela and Cameleopardis 
regions). Since the two reddening maps are independent, 
this argues that the primary source of these residuals is 
due to photometric errors. 

Therefore, while the residual dipole in the SFD compar- 
ison may be real, we conclude that it is more likely due to 
photometric errors. An independent sample is required to 
test this residual dipole more thoroughly. 



an oil slick covered part of the mirror (B87). 
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Fig. 3. — Early-type galaxy extinction residuals versus cosine from dipole direction. The symbol area is proportional to E(B — V)ms,p- 
Crosses indicate galaxies in the Vela region; open circles for all other galaxies. Top panel: SFD E{B— V). Bottom panel: BH E(B— V). 



3. REDDENING ESTIMATES FROM GLOBULAR CLUSTERS 
AND RR LYRAE STARS 

We can perform an independent test of the reddenings 
using the estimated E(B — V) toward distant globular clus- 
ters (GCs) and RR Lyrae stars. 

We take GC E(B — V) estimates from the compilation 
of Harris (1996) 2 , but exclude clusters with |6| < 10 and 
distance perpendicular to the Galactic Plane, \Z\ < 3 kpc. 
This leaves 50 GCs. We find that, in the mean, the esti- 
mated reddenings of the GCs are lower than those of SFD 
by —0.008 ± 0.004 mag. This paper is not concerned with 
zero-point errors (the monopole term) because errors of 
this type drop out of many applications such as peculiar 
velocity studies. We have therefore simply increased all 
globular cluster E(B — V) by 0.01 mag. The scatter of 
E{B-V) is 0.027 mag. 

We have also use the 86 RR Lyrae stars from the sam- 
ple of Burstein & Heiles (1978). The data are reproduced 
in Table 0. Only the first and last four lines of the data 
are given in the printed version. The full table is acces- 
sible electronically. We find, in the mean, that the RR 
Lyrae sample has E(B — V) lower than the SFD predic- 
tions by —0.016 ± 0.003 mag. We have consequently in- 
creased all RR Lyrae E(B — V) by 0.01 mag. The scatter 
around the SFD values is 0.027 mag, the same as found for 
the GCs and for the residuals from the early-type galaxy 
B — V- Mg2 relation. The residuals for these two samples 
are shown in Fig. ^. 

Both the GC and RR Lyrae samples are sufficiently 
small that they do not place strong constraints on the 
residual dipole. Taken together, the best fit dipole is 
F = 0.08 ±0.07 towards I = 312°, b = -47°, although this 
result is not statistically significant; the GC/RR Lyrae 
sample is consistent with no residual dipole. Note, how- 
ever, that this GC/RR F is in the opposite sense to the 
SFD F found above. We can thus use the combined GC 



and RR Lyrae sample as an independent test of the early- 
type B — V- Mg2 residual dipole. We generate 1000 boot- 
strap samples of the globular cluster and RR Lyrae solu- 
tion and project these along I = 80°, b = 23°, the B — V- 
Mg2 dipole direction. We obtain a residual as large as the 
B — V- Mg 2 case only 2% of the time. This argues that 
the B — V- Mg2 dipole found above is more likely to be 
due to a systematic photometric offset between north and 
south. In any case, it establishes the B — V- Mg2 dipole 
result as the upper limit to the large-scale errors in the 
SFD maps. 



4. SUMMARY 

We have compared the SFD and BH reddenings to an 
independent reddening estimate, namely the difference be- 
tween the observed B — V colors of early- type galaxies, and 
the predicted B — V color determined from the B — V- Mg2 
relation. Specifically, we fit a dipole to the residuals in or- 
der to look for large-scale systematic deviations. For BH, 
we find no significant evidence for a dipole once the Vela 
region is excluded. For SFD, we find marginal (92% CL) 
evidence of a dipole in the E(B — V) residuals. We are un- 
able to distinguish between two possibilities: (1) a small 
multiplicative dipolar correction term to the SFD redden- 
ings (a 13% increase towards I = 80°, b = 23°, and a cor- 
responding decrease on the other side of the sky), or (2) a 
small photometric error of 0.01 mag in E(B — V) between 
northern and southern hemispheres. We argue that the 
latter option is more likely, for the following reasons: (1) 
0.01 mag. is well within the expected photometric errors, 
and (2) the dipole for the BH reddenings is very similar to 
that of SFD in both amplitude and direction, even though 
the BH and SFD maps are independent. The sparseness of 
the data do not allow us to determine whether the residu- 
als are actually due to smoothly-varying large-scale (dipo- 
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Table 3 

RR Lyrae and Abt & Golson stars 



I 


b 


E{B-V) 


E(B-V) S FD 


115.70 


-33.10 


0.058 


0.044 


128.40 


-23.60 


0.035 


0.047 


57.90 


-34.00 


0.006 


0.052 


53.20 


-44.30 


0.074 


0.048 


304.60 


57.40 


-0.022 


0.030 


323.40 


67.50 


-0.020 


0.035 


123.90 


18.60 


0.249 


0.247 


122.00 


22.10 


0.112 


0.138 



\ :::■>■ 

.O .■ o • 



Q 



Q °- 



Fig. 4. — Globular cluster and RR Lyrae residuals E(B — V)qc/rr — E(B — V)sfd on the sky in an Aitoff projection of Galactic coordinates. 
Symbols are as in Fig. ^. 



lar) errors in SFD reddenings or in the photometry, or 
whether the problem lies in specific regions, e.g. the Vela 
region and/or the Cameleopardis region. 

The sparser globular cluster and RR Lyrae reddening 
data are consistent with no dipole residual. Nevertheless, 
from this data set, we can exclude the dipole obtained from 



the early-type galaxies at the 98% confidence level. 

Taken together, these results suggest that there are no 
large-scale errors in the SFD maps, or in the BH maps 
(outside the Vela region). The 13% dipole correction to 
SFD should be taken as a reasonable upper limit on any 
large-scale error in those maps. 



REFERENCES 



Burstein, D. & Heiles, C. 1978, ApJ, 225, 40 
Burstein, D. & Heiles, C. 1982, AJ, 87, 1165 

Burstein, D., Davies, R. L., Dressier, A., Faber, S. M., Stone, R. P. 
S., Lynden-Bell, D., Terlevich, R. J., & Wegner, G. 1987, ApJS, 
64, 601 

Davies, R. L., Burstein, D., Dressier, A., Faber, S. M., Lynden-Bell, 
D., Terlevich, R. J., & Wegner, G. 1987, ApJS, 64, 581 

Faber, S. M., Wegner, G., Burstein, D., Davies, R. L., Dressier, A., 
Lynden-Bell, D., & Terlevich, R. J. 1989, ApJS, 69, 763 

Harris, W.E. 1996, AJ, 112, 1487 

Hudson, M.J. 1993, MNRAS, 265, 43 

Hudson M. J. 1999, submitted to ApJ 



Hudson M. J., Lucey J. R., Smith R. J., Schlegel D. J., & Davies 

R. L. 1999, in preparation 
Lynden-Bell, D., Faber, S. M., Burstein, D., Davies, R. L., Dressier, 

A., Terlevich, R. J. & Wegner, G., 1988, ApJ, 326, 19 
Postman, M. & Lauer, T. R. 1995, ApJ, 440, 28 

Press W. H., Flannery B. P., Tuekolsky S., Vetterling W. T., 1986, 
Numerical Recipes: The Art of Scientific Computing (Cambridge 
Univ. Press: Cambridge) 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Smith, R. J., Lucey, J. R., Hudson, M. J., & Steel, J. 1997, MNRAS, 
291, 461 

Trager, S. C, Worthey, G., Faber, S. M., Burstein, D., & Gonzalez, 
J. J. 1998, ApJS, 116, 1 



